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O ! 

^ In the minimal Standard Model (SM) with four generations (the so called SM4) and in "standard" 

two Higgs doublets model (2HDM) setups, e.g., the type II 2HDM with four fermion generations, 
the contribution of the 4th family heavy leptons to the muon magnetic moment is suppressed and 
cannot accommodate the measured ~ 3cr access with respect to the SM prediction. We show that 
in a 2HDM for the 4th generation (the 4G2HDM), which we view as a low energy effective theory 
' for dynamical electroweak symmetry breaking, with one of the Higgs doublets coupling only to the 

4th family leptons and quarks (thus effectively addressing their large masses), the loop exchanges 
of the heavy 4th generation neutrino can account for the measured value of the muon anomalous 
magnetic moment. We also discuss the sensitivity of the lepton flavor violating decays jj, ^ ey and 
i-S^ ' r ^ /17 and of the decay Bs — >■ fJ,fM to the new couplings which control the muon g-2 in our model. 

I. INTRODUCTION 

' Particle magnetic moments provide an important and valuable test of QED and of the Standard Model (SM). In the 

\^ , case of the muon and the electron magnetic moments, both the experimental measurements and the SM predictions 

" are very precisely known. However, due to its larger mass, the muon magnetic moment is considered more sensitive 

[ to massive virtual particles and hence to new physics (NP). 

■ In the SM, the total contributions to the muon 17—2 (aj^*^) can be divided into three parts: the QED, the electroweak 

,—1 ' (EW) and the hadronic contributions. While the QED [1] and EW [2] contributions are well understood, the main 

T-H . theoretical uncertainties lies with the hadronic part which are difficult to control [3| . The hadronic loop contributions 

^ ' cannot be calculated from first principles, so that one relies on a dispersion relation approach Q. At present the 

• ^ , available cr(e"^e~ — )■ hadrons) data are used to calculate the leading<)rder (LO) and higher-order vacuum polarization 

^ ■ contributions to a^*^; the estimated contributions are given by [E,]!] 



PACS numbers; 



'■LO 



6955(40)(7) X 10"", a^fo^"^ = -98(1) x 10"". (1) 



On the other hand, the hadronic light- by- light contribution cannot be calculated from data, hence, its evaluation 
relies on specific models. The latest determination of this term is Q 

afl^" = 116(39) X 10-". (2) 

Including all these corrections, the complete SM prediction is given by 

af*^ = 116591834(2)(41)(26) x 10"" , (3) 
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whereas the current experimentally measured value is [8| 

a;^^P = 116592089(54)(33) x 10"" . (4) 
The SM prediction, therefore, differs from the the experimentally measured value by (see also 

^ncw ^ ^cxp _ ^SM ^ (255 ± 80) • 10~" , (5) 

which allows some room for new physics. For the purpose of this work we are going to assume that the ~ 3ct 
discrepancy in Eq. [5] is due to NP, although we are aware that the estimates of the hadronic contributions have 
appreciable uncertainties that may provide part of the discrepancy. 

In most extensions of the SM, new charged or neutral states^^J, can contribute to the muon anomalous magnetic 
moment (/xAMM) at the one-loop (lowest) level. For example, the /xAMM plays an important role in constraining 
the supersymmetric (SUSY) parameter space, where, as in the SM, the leading SUSY contribution to arises at 
one-loop, and is found to be enhanced for large tan/3. In particular, as was shown in [lol |. SUSY can address the 
observed muon g — 2 discrepancy for tan/3 > 5 and /i > (Higgsino mass parameter), with typical SUSY masses, of 
the particles involved in the loops, in the range 100 GeV — 500 GeV. 

Model independent analysis show that (for details see [9]), for small enough couplings, scalar exchange diagrams 
could account for the observed /zAMM with a scalar mass in the range 480 GeV — 690 GeV, whereas pseudoscalar 
and axial-vector one-loop exchanges contribute with the wrong sign and the one-loop vector exchange contributions 
are too small. 

In this paper we will consider the /xAMM in a new 2IIDM framework with a heavy 4th generation family. Indeed, 
we will show that the ~ 3cr access (with respect to the SM prediction) shown in Eq. [S] can be explained by one- loop 
exchanges of the heavy 4th generation neutrino {v') in a model with two Higgs doublets that we have constructed in 
[ll'] and named the 4G2HDM. These new class of two Higgs doublet models were proposed in [ll| as viable low energy 
effective frameworks for models of 4th generation condensation. In particular, a theory with new heavy fermionic 
states is inevitably cutoff at the near by TeV-scale, where one thus expects some form of strong dynamics and/or 
compositeness to occur. Thus, as was noted already 20 years ago [l2|, the low-energy (i.e., sub-TeV) dynamics of such 
a scenario may be more naturally embedded in multi-Higgs theories, where the new composite scalars are viewed as 
manifestations of the several possible bound states of the fundamental heavy fermions. Besides, our 4G2HDM can 
naturally (albeit effectively) accommodate the large (EW-scale) mass of the heavy 4th generation neutrino, which 
otherwise remains a cause of concern in theories with a 4th family of fermions. 

We recall that an additional fourth generation of fermions cannot be ruled out by any symmetry argument, and 
is not excluded by EW precision data [Ij]. It is also interesting to note, that already the simplest 4th generation 
extension of the SM, the so called SM4, has the potential to address some of the current open questions in particle 
physics, such as the observed baryon asymmetry the Higgs naturalness problem fis'l, the fermion mass hierarchy 
problem [l^ etc....'^' The SM4 can also accommodate the emerging possible hints for new flavor physics [19--24i|. 
However, the SM4 as such cannot explain the observed muon g-2 discrepancy, see e.g., [l^. In fact, even "standard" 
2HDM frameworks (like the type II 2HDM that underlies the minimal SUSY model) with an additional 4th generation 
of heavy fermions, was shown to fail in explaining the measured /iAMM [2^. 

In section|lT]we calculate the /zAMM in the 4G2HDM framework. In sections Hill and ITVl we consider the constraints 
on /xAMM from the lepton flavor violating (LEV) decays /i — >■ 67 and r /i7 and from Bg — > /i/i, respectively, and 
in section |V] we summarize our results. 



II. MUON 2 IN THE 4G2HDM 

At the tree level the muon magnetic moment is predicted by the Dirac equation to be M ~ gfj.^)^S with = 2. 
The effective vertex of a photon with a charged fermion can in general be written as 



u{p')eTf_,u{p) = u{p')e 



2m 



/ 



(6) 



[1] The new states could be a scalar (S), a pseudoscalar (P), a vector (V) or an axial- vector (A). 

[2] Note that the neutral Higgs within the SM4 was recently excluded at the LHC in the range 120 GeV Si mn S 600 GeV [T7| . However, 
this bound is not relevant to a neutral Higgs of an extended Higgs sector, e.g., a 2HDM framework with four generations of fermions, 
such as the ones suggested in [lH . llS| | . 
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FIG. 1: One-loop diagrams for li — >■ Ij^ with charged and neutral scalar exchanges. 

where, to lowest order, Fi{0) — 1 and ^2(0) — 0. While Fi{0) remains unity at all orders due to charge conservation, 
quantum corrections yield ^2(0) ^ 0. Thus, since = 2 (Fi(0) + ^2(0)), it follows that = {g^j, — 2)/2 — ^2(0). 
In our 4G2HDM 11] the one- loop contribution to the /iAMM can be subdivided as 

= K]^^'^^ + (7) 

where [a^Jlp^^^*^ contains the charged and neutral Higgs contributions coming from the one-loop diagrams in Fig.[T] 
(see below; the diagrams with r' and v' in the loop dominate), whereas the SM4-like contribution, [fl/^]^^"*, comes 
from the one-loop diagram with — v' in the loop and is given by [2^ 

[a,\w "^^"^^Aix,,) , (8) 



|C/24p 4V27r 

24 is the 24 element of the CKM-like 

given by 



where t/24 is the 24 element of the CKM-like PMNS leptonic matrix, Xi — mf/M^, and the loop function A{xi) is 



^_ 3a;floga:, , 4a;f - 45a;f + 33x, - 10 
~ (x. - 1)4 + 6(x, - 1)3 • 

For values of m^/ in the range 100 GeV < m^, < 1000 GeV one finds 1.5 x 10"^ < [a^]^^''/|C/24p ^ 3.0 x 10"^, so 
that for |f724p << 1 (as expected) the simple SM4 cannot accommodate the observed discrepancy in a^. 

Let us recapitulate the salient features of the 4G2HDM setups introduced in [ll| . In these models one of the Higgs 
fields (the "heavier" field) couples only to heavy fermionic states, while the second Higgs field (the "lighter" field) 
is responsible for the mass generation of all other (lighter) fermions. Applying this principle to the 4th generation 
leptonic sector we have 



Cy = -El ($^Fe ■(1-1) + ^hYe ■ X) - El [^^lY, -{1-1) + ^hY, ■ Ij vr + h.c. , (10) 

where fL(R) are left (right) -handed fermion fields. El is the left-handed SU{2) lepton doublet and Ye,Y^ are general 
4x4 Yukawa matrices in fiavor space. Also, and are the two Higgs doublets, / is the identity matrix and 
I = diag (0,0,0, 1). The Yukawa texture of (jlOp can be realized in terms of a Z2-symmetry under which the fields 
transform as follows: —^e, ^h, El El, cr -> —cr (for e — e,^, t), ur —vr (for v — , , ) , 

and t'j^^t'j^^v'r^v'r- 

From the point of view of the leptonic sector, the Yukawa interaction in (1101) is the natural underlying setup that 
can effectively accommodate the heavy masses of the 4th generation leptons, by coupling them to the heavy Higgs 
doublet. This setup might also be an effective underlying description of more elaborate constructions in models of 
warped extra dimensions [2^. 

The Yukawa interactions between the physical Higgs bosons and the leptonic states are then given by (see jll|) 

CiMdj) = TT^l \mi^—S^J - • [meXjR + m,^^%L\ \ i,h , (11) 
C{He,e,) - (-m,, + ■ [nii^^R + m,^. Ej^*L] } e,H , (12) 



w 



C{A£,ej) = -ili^l, {m,, tan/3755.j - U ' [miX.R- mi^^*L] } e,A , (13) 

C{H+Uie.j) = Pi { [me^. tan^ • f/,-,; - TOe^//3 • C/feiSL] i? 

v2mvF 
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with 



//3 EE tan/3 + cot/3 , //^ EE — + — , = — - — , 



Sc 



Co 



Sc 



(15) 



and tan /3 is the ratio between the two VEVs. Also, is the charged Higgs, h, H, A are the physical neutral Higgs 
states (h and H are the lighter and heavier CP-even neutral states, respectively, and A is the neutral CP-odd state), 
and £ = e OT ly with weak isospin 1^ = — ^ and = respectively. Also, R{L) = i (1 -|- (— )75) and U is the 4x4 
leptonic CKM-like PMNS matrix. Finally, E'^(E'^) are new mixing matrices in the charged(neutral)-leptonic sectors, 
obtained after diagonalizing the lepton mass matrices 



(16) 



where Lr, Nr are the rotation (unitary) matrices of the right-handed charged and neutral leptons, respectively. Notice 
that and T,'^ depend only on the elements of 4th rows of Lr and Nr, respectively, which we will treat as unknowns, 
i.e., by expressing physical observables in terms of A^fl,4i and Lij,4i or, equivalently in terms of and SJ'^ .t^l 
Following [2^, let us redefine the Higgs Yukawa interactions as 



(17) 



with £ = e or ly and "H ~ , /i, H or A. Then, neglecting terms of order rrie/mr' for e = e, ii,t and terms of order 
Sfj/S^j, for i,j,k = 1,2,3, the above scalar and pseudoscalar couplings, Sj^^,, Sf.^, and Pl^r'^ ^Hu' = or 3), 
which mix the 4th generation leptons with the light leptons, are given in our 4G2HDM by 



4 mw 



pH- 



_g_rn 
2\/2 mw 
9 m 



2V2 m 



w 



tp 



fp^ 



44 



where n° ^h,H or A, = -/^, = ff, F^ = ijp (see Eq. ^ and 



(18) 



TJ* 

TJ* 
(^44 



^4i 
^44 



(19) 



which are the small quantities that parameterize the amount of mixing between the 4th generation leptons and the 
light leptons of the 1st, 2nd and the 3rd generations. In what follows we will take all quantities in Eq. [TS]to be real 
and always set ^44 = E44 = 1 and tan/3 = 1 (for limits on tan/3 in the 4G2HDM see 44^]). We note that and the 
branching ratios for the LEV decays ii — are proportional to 1 -I- cot^ /3 (see Egs. [22land l34p . so that there is no 
enhancement for tan/3 >> 1. 

Using Eq. [THl the charged and neutral Higgs contributions to ag [with — v' and — t' in the loop ("H^ = h,H 
or A), respectively, see diagrams in Fig. [1] are given by (see also |26{ ) 



r 14G2HDA/ 



87r2 



I x{x — 1) 
dx 







oH- 



P"' 



(20) 



r 14G2HDM 



8^2 



dx- 



(1 



p-H" 

flT' 



on" 



pH" 
fir' 



ml, X + m^o(l 



(21) 



[3] Note that since Nn ^i and Lii 4j parameterize mixings among the 4th generation and the lst-3rd generations leptons, we expect 
. < for i, j, k = 1,2,3, see Eq. \W\ 
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Note that, for the neutral Higgs case, the term proportional to m^i /m^ vanishes since \S'^°,\ — \Pl^°i\ (see Eg. [T8)) . 
Therefore, the dominant contribution, by far, to comes from the charged Higgs exchange, in particular, from the 
second term (proportional to m^i iTn^i) in the numerator of Eq. I20( where 



2 



f 13^ 44 



(22) 



so that is proportional to the product Ss^ • ■ 

In Fig. [5] we plot as a function of the product 6-^2 ■ ^U2 (assuming its real) for several values of m^i and 
m}i+ and fixing rrir' = m^i (a^ depends linearly on rrir', see Eg. I22p . Depending on the mass m^/, we find that 
^U2 ' ^^2 ~ 10"'^ — 10^^ is typically required to accommodate the measured value of a^. 




= 500 GeV 



m,.=100 GeV (black) 
=200 GeV (red) 
=400 GeV (blue) 



0.( 



0.016 0.024 0.032 0.04 




H+ = 700GeV 



m,,=100 GeV (black) 
=200 GeV (red) 
=400 GeV (blue) 



0.02 0.04 



FIG. 2: The muon g — 2 as a function of the product 5^2 • Su2, for m^i = 100, 200, 400 GeV, m^i = m^i and with = 500 

GeV (left) and mfi+ = 700 GeV (right). The horizontal lines are the measured 1-a bounds on (see Eq.^. 

In what follows we will consider the constraints from the lepton flavor violating decays £i — £jj and from the decay 
Bs —7- , both of which are sensitive to the quantities du2 and S^^, as will be explained below. 

III. CONSTRAINTS FROM LEPTON FLAVOUR VIOLATION 

LFV decays such as t — > fi'-/ and /i — > 67, which are absent in the SM, are often found useful for constraining NP 
models that can potentially contribute to the fiAMM, as such processes do not suffer from hadronic uncertainties. 
The current experimental 90%CL upper bounds on these LFV decays are [1, [2^ 



Br{T /i7) < 4.4 X 10"* , Br{p 67) < 2.4 x 10"^^ 
Let us define the amplitude for the transition £i — >■ ijj as 

M{e, -> £,7) = Ui^ip') [la^.q'' {A + B75)] u,.(p)e^* , 
where t^* is the photon polarization. The decay width is then given by 



i+!^)(iAp+iBn+4!^(iAp-iBn 



(23) 
(24) 

(25) 



Here again, the new 4G2HDM amplitude M{£i ijj)^^^"^'^''^ can be divided as 



4G2HDM ^ _^SM4( 



(26) 
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where A4^^{£i — )■ ijj) are the SM4-like W-exchange contribution which is obtained from the diagram (right) of 
Fig. [l]with 0^ replaced by plus the diagrams which contain the self-energy corrections to the external fermion 
line £i or ij. In particular, using the definition in Eq. [24] and taking the limit m^^ 0, the net contribution to 
M|^*{e^ £j7) with internal ly' in the loop is given by 29] 



>-w ^ ^0 

aSMA _ nSMA _ eGprnj^ 

Aw - - 4^^2 ^34C/»4-^ {Xu') , (27) 

where Xi — mf/M^ and F(xi) is given by 



_ x^{l~ 6x, + 3xf + 2xf - (3x1 logx,) 
4(l-:r04 ■ 

Here also, we find that — > £^7) is much smaller than the charged and neutral Higgs amplitudes, 
_y^4G2ffDM(-£. ^ ^^j^j _yy^4G2HDM ^ (calculated from the diagrams in Fig.p, for which we obtain 

^AG2HDM ^ ^_[irn,^+rn,^)[sZ:;s^r^,+P^^::P^-^)ir +2rn.\sC ,(29) 



d4:G2HDM 



32t 



A4G2HDM 



pH- 






oH- 


rj- *N 




pU" 




1 



= J,{(-..--..)(5,«rp^>5«>^:;)/r+2m.,(^^^^^ , (32) 

where the loop integrals ^l^^, and /l^' are given by (taking m^. ,m£. <C m|^, m^/ , m^,): 



/f ^ « / dx 



1 ^2 



If"" « j\x 

if « /'dx 

Jo 
/■I 

/2^° K dx 
Jq 



x^{x-l) 





_ a; + m^, (1 - 


^x) 




a;(a; — 1) 






_ a; + m^jjiiX - 


-x) 




a;2(l - a;) 






a; + TO^o(l — 


x) 










a;2 






a; + m^o{i — 


x) 



(33) 

The dominant terms in Eqs. [^M5^ are the ones proportional to m^i from the charged Higgs exchange contribution. 



isZl'P":' - Sf-JP^-^ = -5!l^^2|f,^^|2|5..^|2 L _ (^^^^ _ J ^ (34) 

since the terms proportional to m^' in the neutral Higgs exchanges vanish due to I^J^"/! — jP/'^'l (see Eg. fTB]). 

We thus find that in our 4G2HDM, the decays — >■ 67 and r — >■ /X7 are sensitive to and fea through the 
products (5(72<^Si , <5[/j5s2) and (^[/a&a j '^c/2'^S3)j respectively, so that, in principle, one can avoid constraints on the 
quantities 5ij^ and 5^^ if Jy^, 5ij^^ 5^^ and are sufficiently small. 

In Figs. [3] we plot BR{pL ej) as a function of Sui ■ &2 and ■ fei, for mi,i = 100, 200, 400 GeV, mH+ — 500 
GeV and fixing TOt-' = "t-i^'- We see that for e.g. m^> = 100 GeV and for values of Sjj^ and S^^ of 0(0. 1) [for which the 
product Sjj^ ■ (5^2 reproduces the measured (see Fig. [2])], Sjj-^ and (5sj are required to be smaller than few x 10~^, 
implying that Sui ^ and S^^ ^ 6^^ . 

In Fig. m we plot BR{t — >■ /i7) as a function of • ^s., and (5[/3 • , and in Fig. [5] we give a scatter plot of the 
allowed values in the 6^3 — Sjj^ plane, for which BR{t — > /i7) < 4.4 x 10^^ (i.e., below its 90%CL bound). In both 
plots we use m,' = 100, 200, 400 GeV, mH+ — 500 GeV and we fix nir' — m„' . The individual couplings 5u2 and 
5y,2 are randomly chosen to always be within values that reproduce the measured (see Fig. [5]). We see that the 
products 5u2 ' ^Sa and Sjjz ■ ^2 are required to be at most few x 0.001, in order to be consistent with the current 
bounds on BR{t — ^^). 
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FIG. 3; BR{ii — >■ ej) as a function of Su^ ■ fl^ft) '^'^'^ of ' (fwht), for m^r — 100, 200, 400 GeV, m^i = m^i and 
with mu+ = 500 GeV. Also shown (horizontal line) is the 90%CI/ upper limit on BR{jj, — >■ 67) (see Eg. \23^l . 




0.0015 0.003 0.0045 0.006 0.001 0.002 0.003 



FIG. 4: BR(t — >■ /i7) as a function of 5u^ ■ 5^,2 (left) and 5u2 ■ Ss^ (right), for m^i — 100, 200, 400 GeV, m^i = m,^/ and with 
mff+ — 500 GeV. Also shown (horizontal line) ts the 90%CL upper limit on BR{t — > ^7) (see Eg. \ 23\) . 



We can thus identify a typical benchmark texture for the 4th generation elements of the CKM-like PMNS matrix, 
f7i4, and for the new mixing matrix that can explain the observed /xAMM and still be consistent with the current 
LFV constraints 



U, 



i4 



4i) 



e 

V 1 



(35) 



where e.g., e ~ 0.1 for m^i — 100 GeV. 

Admittedly, the above texture implies a hierarchical pattern which is different from the observed hierarchy in the 
quark's CKM matrix - usually termed as "normal" . Nonetheless, without a fundamental theory of flavor, our insights 
for flavor should be data driven also in the leptonic sector. Besides, the above texture is sensitive to the current 
precision in the measurement of the muon g-2 which can change e.g., if more accurate calculations end up showing 
that part of the hadronic contributions cannot be ignored. 
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FIG. 5: Allowed values in the fes — Sij^ plane, for which BR{t — > fi'y) < 4.4 x 10~* (i.e., below its 90%CL bound), for 
m^i = 100, 200, 400 GeV, m^' = m^i and with mu+ — 500 GeV. The couplings 5u^ and o.'re chosen randomly in the range 
[0,0.2] so that the product Su^ ■ ^2 reproduces the measured a,^ (see Fig.\^. 





SU2 = Ss2 




Br(/i — >• 67) 


Sui ■5t.2 


5u2 ■St.x 


Br(r — >■ /i7) 


5(73. 5e 2 


5v2 -^3 


(GeV) 




xlO^ 


xl0^2 


xlO'^ 


xlO^ 


xlO** 


xlO'' 


xlO* 








1.20 


1.31 


4.57 


1.78 


1.85 


20.40 


100 


0.12 


2.07 


0.64 


2.45 


2.25 


0.77 


3.21 


5.90 








0.14 


1.51 


0.45 


0.26 


2.00 


2.12 








1.21 


0.22 


0.76 


2.03 


0.51 


3.65 


200 


0.05 


2.25 


0.62 


0.51 


0.25 


0.79 


0.41 


2.25 








0.11 


0.43 


0.61 


0.18 


0.16 


1.07 








1.39 


0.16 


0.04 


2.14 


0.73 


0.12 


400 


0.022 


2.19 


0.69 


0.07 


0.10 


0.78 


0.40 


0.21 








0.17 


0.05 


0.02 


0.13 


0.11 


0.14 



TABLE I: The calculated and the branching fractions for the LFV decays r — s- /i7 and /i 67 in the 4G2HDM, for 
several representative values of the couplings 5u2^ ^2 s-nd the products (5(7j,(5e2), (<5i72ifei)i {5u3,5t.2) and (^(721^3) that 
are consistent with the measured /iAMM and which give Br(r — >■ /X7) and Br(/i — >■ 67) at the level of 0(10"^^) and 0(10"^), 
respectively. 



In Table [J we list several representative values of the couplings Su.^, and the products {5ui,5-s:2), {5u2t^t.i), 
{Su3,Ss2) and (Su^tSs^) that are consistent with the measured /lAMM (i.e., a'^P given in eq. H]), and that give LFV 
branching fractions Br(r — )■ try) and Br(/i e-f) at the level of 0(10"^'^) and O(10~^), respectively, that are accessible 
to near future experiments [2^, H^l • 



IV. CONSTRAINTS FROM ^ £+£- 



In the SM, tree level & — ^ s FCNC transitions are forbidden, and also the purely leptonic Bg — >■ i~^£~ decays, 
with £ = e, /Lt, r, suffer from chiral suppression and are therefore very sensitive to new physics. The SM predicted 
branching fractions for these decays are appreciably smaller than those of the semi-leptonic decays. For example, for 
Bs M'*'M~i the SM prediction is [2l| 

Br{Bs fi+n') = (3.2 ± 0.2) • 10"^ . (36) 

In the LHC era the current limit on Br{Bs /i+/i^) has been improved. A combined analysis by LHCb and CMS, 
using 0.34/&^^ and 1.14/6"-^ data sample, respectively, yields [3l| 

Br{Bs fi+i-r) < 1.08 x 10"** , (LHCb + CMS)m5%CL (37) 
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FIG. 6: Dominant SM diagrams in the B^i t decay. 



whereas the same measurement by CDF-II, using a 7/6 ^ data sample, gives [H] 

Br(S^ ^ ^+^-) < 4.0 X 10"*^ CDFmh%CL. 



(38) 



In fact, LHCb has the sensitivity to measure the Br{Bs — > /i^/i^) down to ~ 2 x 10~^, which is about 5cr smaUer 
than the SM prediction. 

In general, the matrix element for the decay Bs (.^(r can be written as [s^ 



M = 



Gpa 



2%/27rsin6'?, 



[Fs U + Fp hd + Fa Pn-i^,-id\ 



(39) 



w 



where is the four momentum of the initial Bg meson and i^^'s are functions of Lorentz invariant quantities. Squaring 
the matrix and summing over the lepton spins, we obtain the branching fraction 



4m^ 



Ml. 



\Fs 



\Fp + 2meFA\ 



(40) 



In the SM, the dominant effect in Bg — )■ £^£ arise from the diagrams shown in Fig. [51 which contribute to Fa in 
Eq. At next-to-leading (NLO) QCD corrections, the net contribution in Fa is given by (s^, IsBj 



I /s. VtbVZ Y{xt) ^-ifp. VtbV;, X 0.997 



mt{mt) 



166GeV 



In the SM4 it is again Fa which receives a non-zero contribution: 

IpSMi / TpSM4\penn . / TpSM4\box 
= {^A ) +[^A ) > 



(41) 



(42) 



where (^pSM4yeng -(-j^g contributions from the Z-penguin diagrams in Fig. [Blwith the top quark replaced by a t' , 
and (^pSM4Yox -(-j^g contribution from the box diagram in Fig. [5] with the replacement t ^ t' and vi v' . At 
leading-order (LO) in QCD, they are given by 



Xf 



i^A ) 



■lfB,VrbV^%\U2i 



6 - Xt: ^ 

1 - Xf (1 - Xt'Y 

^Boixt'.Xy,) 



Inxf 



(43) 
(44) 
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with 



A{1 - Xi){l - Xj) {x^ - Xj){l - XiY 



1 - 2.T, + 

4 



In ; 



{Xj Xi)(l XjY 



1 - 2x, + -4- 



(45) 



In our 4G2HDM, we have additional contributions to Fa coming from the charged Higgs exchange penguin and box 
diagrams (replacing — > H'^ in Fig. [6]) and, in addition, there are new contributions to Fs and Fp. For our 
purpose, we are interested only in the diagrams that are sensitive to the £^i/'H^ vertex, which are, therefore, directly 
related to the muon g-2. The dominant diagrams that contribute to the Bg — ^ £^£~ decay with this vertex are the 
Higgs-exchange box diagrams in Fig. [51 where one or two VF-bosons are replaced by iJ+ and {t, ui) are being replaced 
by both {t,v') and {t' ,v'). Thus, the net contributions to Fs, Fp and Fa in Eq. [39]can be written as 

(46) 



Fs = F^, 



Fp=F, 



H 



and Fa — Fa 



SM 



F 



SM4 



A 1 



where Fg , Fp and F^ are the contributions from the dominant new Higgs-exchange box diagrams, which, for the 
(t' , v') exchange, are given by 



{FsfSl + {Fsf^^H. 



Ml 



nib + nig 



'b 



'"if foH pt 



H- 



H ■ 



oH- 



P, 



H- 



P^7 



pH- 



WH 



{xt',X^>,XH) 



(47) 



F# = (F, 



p)hh 



+ {Fp) 
Ml 



pH- 



qtl r>h 



( pH- oH- 



, P 



BHHiyt',y^' 



{si!;, * + Ff,7 *) } i^t' , , xh) 



(48) 



7H 



box 
HH 



(Fa) 



{Fa) 



box 
WHi 



9' 



{s^fg P^b + Pvs ) {pI^u' S^^jji + 5*^' ^ Bfm{yt' tHu') 



2V2mw iT^w 



-2 \ ^t'bU24 

in/ 



s. 



where Ui = mf/Mfj_ , Xi = m|/M^ and the loop-functions are given by 



BHHi.y^^yi) = Bf^Hivi^yj) 



yi In Vi 



BwHi-'^ii ■'^j T ■'^h) — B^/f^(^Xi,Xj,XH) 



yi-yj l(i-yj)^ (i 

x jj InxH 



yr) 



(i-yO(i-?yj) 



{xH - l)ixH - Xi){xH - Xj) 

x'] InXi 



[Xi ~ l){xi - Xh){Xi - Xj) {xj - l){xj - Xh){Xj - Xi) 



BHHhjuyj) = 



+ 



- Vj)^ {yj - yi) - Vi)^ {yr - yj) 



BwH (Xi ; Xj , Xh ) 



X}j In X}j 



[xH - 1){XH - Xi){xH - X.j) [Xi - l){Xi - Xh){Xi ~ Xj) 



{Xj - l){Xj - Xh){Xj - Xi) 



(49) 



(50) 
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Also, Sj^^, , P/^^, are defined in Eq.[T8]and 



pH- 



{■nib - mt'YiiVt'h + fpVt'b 
{nib + mt')tfjVt'b - fpVt'b 



Vfb 
u* ^tb 
Vt'b 



Vt'b 
Vt'b' 



Vt 



t'b 



(51) 



For the couplings S*// ^ , S"^ , P^, ^ , P^g we use the 4G2HDM Yukawa terms in the quark sector as given in [11[ , 
where E" and E'^ are the corresponding new mixing matrices in the up and down-quark sectors, respectively, obtained 
after diagonalizing the quarks mass matrices. In particular, adopting the type I 4G2HDM of flll|, these matrices are 
given by 



(52) 

in analogy with Eq. I16[ where Dn, Ur are the rotation (unitary) matrices of the right-handed down and up-quarks, 
respectively. They can be approximated by (see jTl| ) 



/O 












\eb\' 



E" 














\ 



(53) 



so that S"''' = if « or j ^ 3,4, and the 34 blocks are parameterized by the quantities and ej. As in [il|, a 
natural choice that we will adopt below is Cb ~ frib/fTib' << £*• Thus, neglecting terms of 0{mb/mt) (and, therefore, 
neglecting also terms proportional to eh), the couphngs S^^ , 5", 



-'t's 



pH- 



P^g can be approximated by 



oH- 

^t'b 



^t'b 



mt' 



Vt 



t'b 



Vt' 



mtfpetVtb 
mtffietVts 



(54) 



leading to Fg — > 



pH 
^ P 



F 



9' 



and 

qH- cH-* 

■^t's '^t'b 
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Ml- 
m^'-mt' 



Re ( S" ,P", 



BHHiyt',y^' 



V2', 



mw 



Notice that the term Un 



{5u2f and {S^^f. 



S^J' - P^.' ) + Vt%U2iS^b {Sf!;, * - P/i: *)}b^j, {xf , X,' , xh)] (55) 

is proportional to (^^2)^ and Re ^^S*^, P^, ^ is proportional to both 
Thus, there is a net effect in Bg — > fj,fi from the charged-Higgs box diagrams even when one 



^Vt'bU24^SP^ 
■ {S^u' 



p": 



of the small quantities that control the muon g-2 vanishes, i.e., when either Sij^ — >■ or 6-^2 
(see Eq.Iini)- For example, for ^ we have [a^]^^^^*^ 0, while 



0, for which cases 



,„ UG2HDM 



{Su2 ^ 0) 



'^t's '^t'b 



-BU{vt'.V^^')^fl\U,A^\T.l,\^\5^2 



(56) 



The contributions from the charged Higgs exchange diagrams with {t, v') can be obtained directly from Eqs. I47I48I 
and 1491 by replacing t' with which eventually requires the replacements: t' ^ t and E^^g in Eg. 1511 

In Fig. [7] we plot the contribution to BR{Bs /i/i) in the 4G2HDM from the box diagrams with the and {t, v'), 
{t', v') exchanges, as a function of et (fixing Cb — 0.01) and A^^ = VfbV^^, for mi,' = 100, 200, 400 GeV, m,-/ — m^/ 
and mfj+ = 500 GeV. The allowed ranges of the key parameters 5u2 and 8^,2 , which control the muon 5 — 2 in our 
model, are randomly chosen in the range [0, 0.2] to be consistent with a'jf^ given in Eq. S] We also show the current 
experimental bound [3l| and the SM predicted value for BR{Bs ^/i). We see that the contribution from the new 
(i.e., in the 4G2HDM) box diagrams that involve the heavy 4th generation neutrino is consistent with the current 
experimental bound on BR{Bg — >■ fifi) for values of (5[/2 and 6-^2 that reproduce the observed muon g — 2. It is also 
interesting to note that both in the SM4 and in the 4G2HDM, Bt{Bs — )■ /i+/i~) can differ from the SM value by 
at-most a factor of 0(3) in either direction. 
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FIG. 7: BR{Bs j-ij-i) m the 4G2HDM from box diagrams with the and {t, i/'), {t' , v') exchanges, as a function of et (fixing 
€b = 0.01^ and = Vt'bV^i^, for ra^i — 100, 200, 400 GeV, m^i = m^i and Tnu+ = 500 GeV. We considered only values 
of '^^d, 5y,2 which are allowed by given m Eq.\^ keeping both of them ^0.2. Also shown are the experimental 95% CL 
upper bound (upper /red horizontal line) and the SM predicted range (la) of values (lower /black horizontal lines). 



V. SUMMARY AND DISCUSSION 



We have considered the effects of 1-loop exchanges of heavy 4th generation leptons on the muon g — 2, on the 
lepton flavor violating decays fi ej, t and on Bg — ^ in the 4G2HDM which is a 2HDM where the 

Higgs doublet with the heavier VEV is coupled only to the 4th generation doublet while the "lighter" Higgs doublet 
is coupled to fermions of the lst-3rd generations. This model is particularly motivated for the leptonic sector, as it 
effectively addresses the heaviness of a 4th generation EW-scale neutrino. 

The muon g — 2 is sensitive in our model to the product ■ ^y,2 , where Su2 = , ^t,2 = , Uu is the leptonic 
CKM-like PMNS matrix and E?,-, are new mixing matrices in the charged and neutral leptonic sectors that are 
unique to the 4G2HDM. 

We find that, depending on the mass myi, the experimentally measured muon magnetic moment can be accounted 
for if 0(10^'^) £ 5u2 • ^2 ~ O(10~^). We also find that the decays ^ — >■ 67 and r — 7> ^7 can have branching ratios 
which are not too far below the current bounds, i.e., of 0(few • 10~^^) and ©(few • 10~^), respectively, if the products 
5ui ■ fea: ^U2 ■ fei ^ C(10~^) and 5u2 ■ ^Ss, ^u^ ■ ^2 ^ 0(10"^), respectively. 

We also considered the effects of one-loop exchanges of the 4th generation heavy neutrino v' on the decay Bg M^A*^ 
and found that, in the four generations model considered here, BR(i3s — > can be larger or smaller than the 

SM predicted value by a factor of about three, for values of 0(10"^) ^ Su2 ' ^2 ^ O(10~^), which render the observed 
value of the muon g — 2 to be consistent with the current upper limit on this decay. 
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